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ABSTRACT 


Calculating the lift, drag, and pitching moment coeffi- 
cients for an airfoil from the static pressure distribution 
obtained from wind tunnel tests is routine task when steady 
flow is considered, but it is much more complicated when the 
airfoil is operating in an unsteady flow field, similar to 
that experienced by a helicopter rotor blade, produced by 
an oscillating wind tunnel. A data reduction routine capable 
of condensing the large numbers of data associated with the 
unsteady investigation, as well as a numerical integration 
algorithm for the unsteady aerodynamic coefficients, were 
developed; however, no unsteady data were collected due to 
hardware failures. The ability of the program was demonstrated 
on previously obtained steady and quasi-steady data and sampie 


Be ults were presented. 
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I. INTRODUCTION 


The primary objectives of this investigation were: 
1) to report on the data reduction techniques used to ob- 
tain values for both steady and non-steady aerodynamic 
coefficients, resulting from the experimentally acquired 
static pressure distribution over the model when the tunnel 
was operating at a specified dynamic pressure, with Coanda 
blowing oscillated at various frequencies and amplitudes of 
momentum blowing coefficient; and 2) to present a sample 
of results. Removal of the model from the tunnel for 
modifications in November 1976 precluded any actual process- 
ing of non-steady data consistent with the stated goals.” 

However, the data reduction techniques are applicable 
to the pressure distribution obtained with steady Coanda 
blowing and will be demonstrated on such data obtained in 
April 1976. Prior to removal of the model, the tunnel was 
operated without blowing in an effort to generate a set of 
Mon steady data to be used in tailoring the data reduction 
procedures to be used downstream in the on-going project. 
In that test, the model experienced self-excited harmonic 


pressures resulting from vortices being shed alternately 


Ihe modifications, to correct for surface delamina- 
tions near the upper jet nozzle, include a steel upper 
surface trailing edge and slot height adjustment to 0.015 
inches. Additional pressure taps were installed on the 
upper surface forward of the Coanda slot, and leaking 
tubes were repaired while the modifications were being 
implemented. 





from the upper and lower surfaces. The manually recorded 
data consisted of the RMS and mean values of the unsteady 
pressure, and the phase angle in respect to a reference for 
each pressure orifice. A more complete description of the 


data acquisition system is given in Section II B. 
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Il. EQUIPMENT AND INSTRUMENTATION 


A. AIRFOIL 

The airfoil used in this study was obtained from the 
Lockheed-California Company after completion of their 
design feasibility study on circulation control rotors 
(Ref.1). It had an elliptical cross-section with a chord 
of 10.215 inches, thickness ratio of 0.20, and camber of 
3.3% (Figure 1). A slot to facilitate Coanda blowing was 
machined in the upper surface at chord position X/C=0.95, 
with the slot height set at 0.010 inches. Coanda air was 
supplied to the plenum cavity at both ends of the section 
to insure uniform flow along the span. Fifty-three 
pressure taps at mid-span with a chord-wise orientation were 
monitored through two scanivalve transducers for an elec- 
trical representation of the pressure distribution over the 
airfoil. The pressure tap locations are listed in Table I. 

The 2.0 x 2.0 foot unsteady flow wind tunnel in the 
Department of Aeronautics Laboratories at the Naval Post- 
graduate School was used as it is capable of simulating the 
velocity flow field experienced by blade elements of a 
helicopter in forward flight. The model was mounted at an 
angle of attack of -5 degrees, which closely corresponds to 


the zero-lift angle of attack for the data reported herein. 


ПЕД, 
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AIRFOIL PRESSURE THP LOCATIONS 
ТАР + 
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В. INSTRUMENTATION 

The data acquisition system evolved from a classic 
manometer board arrangement with readings hand recorded 
(Ref. 2) and reduced to the present system, which employed 
electrical read-outs with an existing Digitec system with 
customized signal conditioning. The digital multiplexing 
circuitry in the Digitec system was interconnected through 
the scanivalve solenoid controllers for automatic scani- 
valve advancing as the output of the digital voltmeter was 
sequentially printed on a paper tape. Pressure sensing was 
accomplished by connecting each of the pressure taps to 
one of the two scanivalves (Figure 2) through tubing of 
uniform diameter and length. This was done to insure a 
constant transfer function for all pressure taps, when 
measuring unsteady pressures. Calibration tests using the 
method reported by Johnson (Ref.3) showed an acceptable 
variation in dynamic gain in the frequency range from 0 to 
300 Hz (Figure 3), which included the range of practical 
interest. 

For the unsteady aerodynamics analysis a the 
voltage representation of the surface static pressures was 
of the form 

e =a tb (t). (2:1) 
Full analog measuring techniques were employed in measuring 
the mean voltage, RMS voltage, and phase relationship to a 
reference for each pressure tap. The mean voltage was ob- 
tained by processing the signal through a low-pass filter 


circuit (Figure 4) with a two second time constant which 
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Airfoil Instrumentation 
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earresponds.te a comer frequency of 0.08 Hz on a Bode plot. 
Therefore, the influence of the unsteady voltage was essen- 
tially negated and the output readings represented the 


steady components of the signal. 


C. ADVANCED DATA ACQUISITION SYSTEM 

To meet the primary objective mentioned in Section I, 
a more sophisticated data acquisition and reduction system 
will be required and is now being constructed by Englehardt 
at the Naval Postgraduate School in the Department of Aero- 
nautics. The system, based on the Intel MDS-80 micro- 
processor unit, will control the logic of the experiment, 
collect and store data, and handle the entire data reduction 
process. The algorithm for data reduction discussed in 
Section V was proven sufficiently simple and capable of 


performing all required data reduction and numerical integra- 


tion. 


“Master of Science thesis by Englehardt with a complete 
description of the entire microprocessor system will be 
published in June 1977. 
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III. DATA REDUCTION 


A. STEADY DATA PROCESSING 

The aerodynamic coefficients are defined by integral 
equations but must be approximated by numerical integra- 
tions, since data are available only at a finite number 
of pressure tap locations. The approximate steady normal 


force, chord force, and pitching moment coefficients are 


1.0 
с = | (Ср. -Ср 24070) 
Y/C (max) 
с. = ATI) (5.1) 
ао) 
Y/C (max) 
— (6, IE 
-Y/C (min) 
0 x 5 
| Cy а-а, 


Then, with an angle of attack correction for the force 
coefficients and a moment transfer to the quarter-chord 
position, the aerodynamic coefficients are 
C,=C cosa-C sina 
= 1 - 
Са C, Sina Ce cosa (34,2) 


Enle/s) “Сп(тЕ) 27264. 


The pressure data from the April 1976 tests and the 


November 1976 mean pressure data were transferred to a 


NS 





Hewlett-Packard 9830 tape cassette and processed using the 
program listed in Appendix A, which employs the equivalent 


numerical relationship to the equations stated above. 


В. UNSTEADY DATA PROCESSING 

Before the unsteady data were suitable for numerical 
integration, a dynamic gain correction was necessary to 
account for the transfer function discussed in Section II. 
The frequency of the self-induced oscillations was constant 
at 262 Hz when the tunnel airspeed was approximately 100 
ft/s. From Figure 3, the corresponding gain correction, 
Gel, was approximately 0.5, which implied that the amplitude 
of the surface static pressure was attenuated by a factor 
of one-half at the transducer. All phase angle values were 
relative to tube number 7, which was arbitrarily chosen as 
a clock reference. Phase angle errors resulting from the 
two different scanivalve circuits were faired out in the 
analysis (Figure 5). 

With the above corrections, the numerical integration 
of the unsteady static pressure distribution was carried 
out in a manner similar to the steady pressure integrations. 


For the unsteady normal force coefficient 


1.0 
3 
9 sinó = | (с, sine-C, 81092407 
0 
1, 1.0 
(С > 059 = | (С_ со<ф-С совф)а(Х/ ) 
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Then “2 = (€, ^ sine) ^«( €, ^» cose)” 


227) 
And фт = tan”? СС. ^ sine) 
een 
L 
( Es cosó) 
1 


3 
where i 15 the RMS value ot the normal force coefficient 
and j+ is the phase angle relative to tube number 7. The 


chord force and pitching moment coefficients were calculated 


similarly. 
y Y/C(max) 
er sing = (C, suno t, sine)d( /Q) 
-Y/C(min) t Е 
1 Y/C(max) 
2 ү 
(c ^) cose = (С, <о59-С, соз%)4 7/0) 
-Y/C(min) Í 5 


1, 
1 É : шонг 
E^ = (€ °) sing) +(<C “5 con : $T = E > $1Пф) 


(4,7 y cose) 


1, Y/C (max) 
C 2)” sing = (С. віпҙ-С біп) С/оас/оз 
S f Pr 
-Y/C(min) 
1.0 
| <» sing-0, sino aa 
0 u 
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The unsteady aerodynamic coefficients are calculated as in 
Eqs. 3.2, with the aid of elementary trigonometric identities 
account for the phase angles. The program listed in 
Appendix A has the option for unsteady numerical integration 
incorporated. 

i IS intended that the pap oscillating pressure 
measurements be repeated when the improved data acquisition 
system becomes available. The tests will be conducted at 
several tunnel dynamic pressures to determine the Reynolds 
number dependence of the data. The previous measurements 


were at approximately R,70.54 X S. 
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lV. FT RESULIS 


ПІ STEADY DATA RESULTS 

The validity of the concept of circulation control was 
demonstrated by Englar (Refs. 4 and 5) and is confirmed 
here with the steady momentum blowing coefficient results. 
The resulting aerodynamic coefficients are summarized in 
Table II. The ability to produce relatively large in- 
creases in lift coefficient with only a small increase in 
blowing coefficient, shown graphically in Figures 6a, b, 
and c, is one of the desirable characteristics of the 
system. Probably the best measure of the effectiveness is 
the aerodynamic lift amplification ratio, 907/30, which 
was computed to be approximately 15.3, the low value being 
attributed to the geometric deficiencies of the trailing 
edge slot as alluded to in Section I. However, the results 
are encouraging. 

Close scrutiny of the coefficient of pressure informa- 
tion plotted against X/C (Figures 7-12) reveals a pressure 
distribution over the model such that the center of pressure 
was at the mid-chord position, which is in-agreement with 
the theory. The center of pressure position is also 


Verified by the graph of Cn(c/4) Vs. С 


u 
Ba CS Las 
CP. (30,73€. ) (15.5) 
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TABLE II 
AERODYNAMIC COEFFICIENT SUMMARY 
STEADY AERODYNAMIC COEFFICIENTS 


Angle of Attack--5 degrees U=100 ft/s 


en 





UNSTEADY AERODYNAMIC COEFFICIENTS 


С. --0.0765 C4 (RMS) 22.1404 $2225 degrees 
C¿= 0.0451 C,(RMS)=0.1629 ф= 39 degrees 
C, "70-1065 C, (RMS) =0. 5817 b= 63 degrees 


25 





1.00 
.Вӣ 
‚ER 


Ú 
B 
В. ЧИ 
B.22 
D 


‚Ид 


4.08 
Й.ЙБ 
B.2H 
4.42 


Й.ИЙ 


=Й ИЙ 


= 
ы 


В. ИБ 


FIGURE ЕН 


Mm 
E 
tJ 


И. ИБ 


FISURE БЕ 


m 
гы 


се 
FIGURE BC 


Й. ИБ. 


STEADY AERODYNAMIC COEFFICIENT 
VARIATION WITH Co 


26 


1.23 








ZDEFFICIENT OF PRESSURE 


(P Vs. X/< 


заза 
RUN ND. 42341 
=Й. ИИ 
-Ч.й 
UPPER SURFACE 
-3./ 
-2.2 
E. 
A + + + + + EE 
+ + 
+ ? + 
0.2 | Х/С 
+ 22 d. Й.Б Й.В 1.4 
+ 
ED 
+] . Ú 
++ 





Figure 7 


211 





«ЛЕ ЕЕГ СЕЧИ ЕЕ PRESSE 


Ld 


cp VS. WC - 
RUN ND. 42322 


和 = 有 .BBSS 


UPPER SURFACE 


4.2 1.1 Й.Б B.B 





28 





<IEFFI<Z TEN T ae PRESSURE 


Сы руда 


RUN ND. 42343 
=Й.ИБЧИ 


UPPER SURFACE 





риге 9 


= 





ӘБПЕЕКТСТЕРТГ OF егес ЕЕЕ 


= 


= 


(PNS А/С 


RUN NO. 42344 
(и=Й.ИЧЧВ 


UPPER SURFACE 


: + + ++ Ж : 
+ 
+ 
+ 
+ + 
#Х/С 
4.2 日 .4 Й.Б B.B L.A 





Figure 10 


50 





COEFFICIENT OF PRESSURE 


- 3 


CP V5. XC 
RUN ND. 42385 


(и=й.И272 


UPPER SURFACE 


2.2 d. B.E 2.8 





Figure 11 


Sl 





COEFFICIENT OF PRESSURE 


(P V5. XC 


-5. 
RUN NO. 42906 
C 24.4133 
-W.p 
UPPER SURFACE 
-3.2 
-2.0 
+ 
-1.р И + 
А NEED COMER. 
+ + 
+ + 
B. i 7 
. Й.2 g.u Й.Б J.B 2 
+ 
| 
+1.2 





Figure 12 


52 





Also, the low drag qualities of the airfoil are favorable 
when full thrust recovery of the momentum blowing coefficient 
is assumed. Pressure coefficient data are listed in 


Appendix B. 


В. UNSTEADY DATA RESULTS 
The unsteady aerodynamics of the elliptically shaped 
airfoil have been given a rather cavalier treatment in 
the past but their contribution to the total forces and 
moments must be realized. The implications of the unsteady 
aerodynamics is not generally well understood; therefore, 
a simple example is presented here to emphasize their 
importance before discussion of the unsteady results. 
Consider an airfoil at zero angle of attack with the 
RMS pressure distribution equal on the upper and lower 
surfaces (Figure 13a) but with unequal relative phase 
angles (Figure 155) taken from the quarter-chord and the 
three-quarter-chord positions from the actual unsteady data 
(Figure 5). The unsteady normal force coefficient, с 2), 
is calculated from the numerical approximation to Eq. 3.3, 
as: 


% 
E^ sinó$- 0.5(sin(210°)+sin(230°)-sin(35°)-sin(55°)) 


25 
(c, ^) sine = -1.33 


0.5 (соѕ (2109) +соѕ (230°) -соѕ (359) -соѕ (55°)) 


1. 
E ) cose 


L 
аы ) coss = -1.45 


(<) = 41.332 + 1.45 б = Te 
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(2) = 1.97 Ше 
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Similarly, for the moment coefficient (referenced to the 


trailing edge): 
i 
€ 2) sine = 0.5((sin(210°)-sin(35°))0.75+(sin(230°) -sin(55°)).25) 
1, 
Du a 
С) $1пф = -0.6 
2 É > O O O O 
(C DEC. oe соз 2105) -со5(25 ))0:- 75ғ(со05(250”)-со5(55 ))0.25) 


2. 4 
€, ) cos = -0.78 


1 
7 CE 2 2 _ -0.6 
(C, j"- Jor6^ - 0.78 $ = 0 75 
2 
Ы 2 О 
оо ф 2117 


Clearly, the force and moment coefficients are of significant 
magnitude even if a uniform distribution (RMS) of pressure 
coefficient were assumed, because of the importance of con- 
sidering relative phase relationships. The dynamic behavior 
of the unsteady aerodynamics of the circulation control 
airfoil, particularly the lift amplification ratio, in the 
unsteady environment caused by cyclic blowing will be in- 
vestigated after the improved data acquisition system is 
operational. | 

As discussed in Section I, measurements were taken in 
November 1976, with the airfoil experiencing harmonic 
pressure variations. The unsteady data are listed in 
Appendix C. The experimental data were numerically inte- 
grated, similar to the previous example, for the lift, drag, 
and pitching moment coefficients (Figures 14 and 15). It 


should be recognized that the phase shift between the upper 
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and lower surfaces is approximately 180 degrees, the result 
being the relatively large amplitude of the aerodynamic 
coefficients. Highlights of the unsteady results, summarized 
in Table II, were that: 

l. The time varying pressures over the model were 
predominantly harmonic, with a non-dimensional frequency 
(Strouhal no.) of 


E Е 
S = =~ = 0.449. 


2. The RMS value of the normal force coefficient was 
approximately 2.14, which was an order of magnitude greater 
than the RMS drag coefficient. 

3. The test Reynolds no. was approximately 


R, = 0.54 x 10°. 


4. The location on the airfoil for the center of 
pressure due to unsteady loadings was approximately 0.491 
chord and the RMS value of pitching moment at this reference 
point was 


1 
2 


(C4) = 0.1694 
m 


with a phase angle of 90 degrees lag relative to the normal 
force. This term may be considered as an apparent mass 


pe effect. 
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V. DATA REDUCTION FOR THE ADVANCED DATA REDUCTION SYSTEM 


The advanced data acquisition system described in 
Section II C will allow a more comprehensive investigation 
of the unsteady aerodynamics of the model, but the data 
reduction techniques must be tailored to handle the large 
numbers of data to be collected. Oscillation of the plenum 
pressure as described by Bauman (Reference 6) induces 
harmonic blowing from the trailing edge slot at a specified 


frequency, W It is assumed that the unsteady static 


17 
pressure at each station will have a steady value plus a 
periodic value with frequency Wi superimposed due to the 
oscillating momentum blowing coefficient. 

Since the data are assumed periodic in form, the un- 


steady pressure at any point on the airfoil may be defined 


by a Fourier series (Ref. 7) as follows: 
Р(х,%)-А.(х)% 24 [А (х) cos (nW,t)+B_ (x) sin(nW,t) (2521) 


where W. is the frequency of the oscillation. The standard 
formulae for obtaining the Fourier coefficients for the 


infinite series are: 


T 
AQ) =e | P(x) dx 
Zn 
T 
N = + | P (x)cos (nx) dx (2522) 
-1 
1 (4 | 
B (x) = = | P(x)sin(nx) dx 


E 





Because the experimental data are in the form of a 
discritization of a truncated periodic signal, the above 
equations are slightly altered for estimating the Fourier 


coefficients. The modified form is: 


K 
À EE 
“оог ES ; (x) 
x 2; 
а) = ү E. P;(x)cos(nW,t) ‚ЭО, уп (5,5) 
% 2 К 
Zu) = ү га P;(x)sin(nW,t) 


where k is' the number of samples to be analyzed and n 
15 the harmonic number. It is recognized that the physical 
meaning of A (x) is the steady or mean value of the oscillat- 
ing pressure. That is, A. (x) = (py G2) where 
| Т 
(PL) = 12" 3 | P(x) dx (5.4) 
-T 
By subtracting out the mean pressure value, the signal becomes 
biased to a pure oscillating pressure with a mean value equal 
to zero. The terms Á (x) and В (х) are equivalent to the 
cosine and sine components, respectively, discussed in 
section III B. Figure 16 illustrates the reconstruction of 
a discritized signal utilizing this algorithm programmed on 
the Hewlett-Packard 9830 calculator (program listing in 
Appendix A). The error was minimized by insuring that the 
number of samples analyzed was an integer multiple of the 
frequency; however, the truncation error was never greater 


than 23. 
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Having developed an effective method for handling the 
large numbers of data, calculation of the aerodynamic force 
and moment coefficients can be accomplished. Before the 
actual processing of data begins, the sample pressures are 
converted to non-dimensional coefficient form to expedite 
calculations. The unsteady pressure coefficients are 

NO = i= 1,2,3,...,53 (5.5) 
where P and H are steady pressure values from a pitotstatic 
tube mounted on the floor of the wind tunnel. Then the 
oscillating signal at any pressure tap is; 


C (x. t) = 2 ¿e cos (nw, t)+B_ sin(nw,t) 


E 05 = у T, sin(nw,t+9_) (5.6) 


where Т = NA в 1 > 9 - “ву, 

n 
The unsteady тэ are related by the phase angles y , 
however, the phase angles are meaningless unless they are 
measured relative to a "clock" or reference signal. There- 
fore, it 1S necessary to designate one pressure tap to be 
used as the reference. Then, the relative phase angle is 
227779) 


ф = Q. - % 


rel П ref 
TOT simplicity, $; shall hereafter be interpreted as the 
relative phase angle. 

With the steady pressure coefficient removed from the 
oscillating signal, separate aerodynamic coefficients can 


be calculated for the steady and unsteady signals. The 
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steady force and moment coefficients are found exactly as 
in Section III A. The unsteady aerodynamic coefficients can 


be defined in a Fourier series form as: 


C(t) = E “ny sin (kw, t+, ) 
C (t) = y e sin (kw, t+9,) (558) 
C CO = E EN Sin (kw, t+, ) 


The magnitudes and phase angles of the coefficients can be 
determined thusly: 
for the first harmonic 


E m = В 


where 1 
m e sino. (5-c (sim (Sha 
p" - | (с, (езіле) E (sins H) de) 
0 u 
С соѕф = i E (Х) сов (>) - C (5) cos (2) d (X) (579) 
n; ? py `C ШІ С Doc Pu >) С i 
0 u 
then 
_ C_ sing 
С = (6 cos+,) + (C sin#,) , $1 - Tan 1 id 1 
121 ai x Cn С059, 


1 
Similar calculations.are done to obtain values for the higher 
Ec. ) as well as the chord force and 


1 
pitching moment coefficients using the higher harmonics of 


harmonics Cen бал? 


pressure coefficient information- as obtained during the 


initial data reduction previously described. 
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VI. CONCLUSIONS 


As mentioned in Section I, the two-fold objective of 
this paper was to develop a computer program for unsteady 
aerodynamic analysis and to demonstrate its capabilities 
on experimentally generated unsteady data. The objective 
was partially satisfied; a computer program, well within 
the limitations of mini-calculators such as the Hewlett- 
Packard 9830, was developed which employed a straight- 
forward numerical integration of pressure coefficients to 
obtain values for the aerodynamic coefficients. The pro- 


gram can be modified for unsteady data by properly taking 


into account the cosine and sine components of the pressure 


signal relative to a fixed harmonic (clock) signal. Also, 


the extraction of the Fourier coefficients of the truncated, 


time-varying signal appears to be within the capabilities 


of the updated data acquisition system. 
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APPENDIX A 
COMPUTER PROGRAM LISTINGS AND OUTPUT 


REM PROGRAM FOR DECOMPOSITION OF A FUNCTION 
REM INTO FOURIER COEFFICIENTS 

DIM EC256] 

W1=(20*PI) 

DISP "DIVISIONS PER CYCLE?"; 

INPUT D 

P=2*PI/W1 

T=P/D 

REM COMPUTE SAMPLE FUNCTION 

FOR I=1 TO 220 

X=W1* (1-1) *T 
ECIJ=SIN(X)+0.5*SIN(2*X)+0.3*SIN(3*X) 
NEXT 1 

PRINT 

PRINT "THE DRIVING FREQUENCY IS 10 HZ" 


PAINT D; “DIVISIONS PER CYCLE" 
PRINT X 

PRINT " COEFFICIENTS" 

PRINT 


PRINT ТАВ2"С05"ТАВ15"51М" 

FOR J=1 TO 3 

PRINT 

PRINT TAB2"A";J, "В",7, "PHASE ANGLE" 
А1-81-0 

FOR I=1 TO 220 
Bl=(B1+E[LIJ*SIN(J*W1*(I-1)*T)) 
Al=(A1+ECI] *COS(J*W1* (I-1)*T)) 

NEXT I 

P=ATN(B1/A1) 


WRITE (15,170) (A1*2/220), (B1*2/220), (P*180/PI) 
FORMAT F6.3,7X,F6.3,13X,F5.1 
NEXT J 
END 
THE DRIVING FREQUENCY IS 10 HZ 
10 DIVISIONS PER CYCLE 
COEFFICIENTS 
COS SIN 
A 1 Baal PHASE ANGLE 
0.000 1.000 90.0 
A 2 B 2 PHASE ANGLE 
0.000 0.500 90.0 
A 3 B 3 PHASE ANGLE 
0.000 0.300 90.0 
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130 
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150 
160 
170 
180 
190 
200 
201 
202 
#05 
204 
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206 
210 
2117 
213 
214 
216 
2117 
218 
220 
222 
224 
226 
227 
228 
230 
255 
240 
242 
244 
245 
246 
247 


REM COMPUTE CL,CD,CM--STEADY AND UNSTEADY 
REM FROM CP DATA ON TAPE CASSETTE 
COM R1,D[80] ,P[80] ,C[60] 

DIM AC106] 

PRINT 

LOAD DATA 4,A 

С9-0 

Х9-0.95 

Y9=0.048 

DISP ''CP DATA FILE NO. IS"; 

INPUT Z 

LOAD DATA Z 

рС221= 0.4 

рС361=0.539 

61=1 

PRINT "INPUT OPTION" 

PRINT "1 FOR STEADY DATA" 

PRINT "2 FOR NON-STEADY DATA" 

INPUT G3 

GOTO G3 OF 130,212 

GOTO 130 

PRINT 

PRINT ''RUN NO. IS"R1 

PRINT TAB7"UPPER"TAB24"LOWER" 

PRINT TAB2"N'TABO''CP'TAB19"N'TAB26''CP" 
FORMAT F3.0,2X,F7.3,5X,F3.0,2X,F7.3 
FOR N=1 TO 27 

M=N+26 

WRITE (15,170)N,C[N] ,M,C[M] 

NEXT N 

FIXED 3 

EM M "UGPIS"UCO"AT X/Cs"XO". Y/C="Y9 
STANDARD 

DISP "CORRECT THE CP VALUES" 

STOP 

GOTO 270 

DISP "INPUT DYNAMIC GAIN"; 

INPUT G2 

FIXED 3 | 

PRINT "RUN NO. IS'R1 

C9=C [21] +0.33* (C [22] -C[21] ) 

PRINT "TAP"TAB7"CP(MEAN)"TAB17"CP(RMS)"TAB27"PHI" 
FORMAT F3.0,3X,F7.3,3X,F7.3,3X,F4.0 
FOR I=1 TO 53 

WRITE (15,220)1,C[1),D[1J],P[1) 

NEXT 1 


PRINT "C9="C9 

MSP" CORRECT THE CP VALUES" 
STOP 

GOTO 270 


FOR I=1 TO 53 

C[I]=G2*D[IJ *SIN(P[I1/57.296) 
NEXT I 
C9=(D[21]+0.33*(D[22]-D[21])) 
G1=2 

GOTO 270 
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250 
252 
254 
256 
257 
258 
270 
272 
274 
276 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 


FOR I=1 TO 53 
СГІ1-С62%рГІ1%С05(РГ11/57.296) 
NEXT I 
51212211 50:-532:(81221-0(211)) 
61-3 

GOTO 270 

L1=0 

L2=0 

Р1-0 

Р2-0 

FOR I=1 TO 26 
IF I=21 THEN 490 
C1=-C[I] 
C2=-C[I+1] 
X1=A [1I] 

X2=A[ 1+1] 
GOSUB 1120 
L1=L1+N1 
P1=P1+M1 
C1=C[I+27 
C2=C [1+26] 
X1=A(1+2 7] 
Х2-АГІ%26) 
GOSUB 1120 
L2=L2+N1 
P2=P2+M1 

NEXT I 

GOTO 640 
C1=-C [21] 
C2=-C9 
Х1-А(211 
X2=X9 

GOSUB 1120 
L1=L1+N1 
P1=P1+M1 
C1=-C9 

E2=-C| 22] 
Х1-Х9 

Х2=А 22 

GOSUB 1120 
L1=L1+N1 
P1=P1+M1 

GOTO 400 

D1=0 

D2=0 

P3=0 

P4=0 

FOR I=1 TO 13 
221-0(11 
С2-0(Т1-11 
Ү1-АГІ»521 
Y2=A(I+54] 
GOSUB 1160 
D1=D1+N1 
P3=P3+M1 

NEXT I 
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770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
950 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1050 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
#120 
150 
1140 
L150 
1160 
1170 
1180 
IIS 
1200 
1205 
12109 
152209 
229 
1230 
1250 
1260 
1270 
1280 
1290 
1500 


FOR I=40 TO 52 

ОШО | 

C2=C[I+1] 

Ү1=А[1+53)] 

Y2=ALI+54] 

GOSUB 1160 

D1=D1+N1 

P3=P3+M1 

NEXT I 

FOR I=14 TO 39 

IF I=21 THEN 970 

C1=-C[1+1] 

созе 

Y1=A[LI+54] 

Y2=A[1+53] 

GOSUB 1160 

D2=D2+N1 

P4=P4+M1 

NEXT I 

GOTO 1200 

С1=-С9 

C2=C 21 

Ү1-Ү9 

Y2=A 74 

GOSUB 1160 

D2=D2+N1 

P4=P4+M1 

С1=-СЁ22] 

С2=-С9 

Y1=A[ 75] 

I Y 

GOSUB 1160 

D2=D2+N1 

P4=P4+M1 

GOTO 950 
M1=C1*(X2-X1)*(1-0.5*(X2+X1) 
M1=M1+(C2-C1)*(X2-X1)*(0.5-(X2/3)-(X1/6)) 
N1=0.5*(C2+C1)*(X2-X1) | 
RETURN 
N1=0.5*(C2+C1)*(Y2-Y1) 
M1=0.5*C1*(Y2-Y1)*(Y2+Y1) 
М1-М1%(С2-С1)%(Ү2-Ү1)%((Ү1/6)%(Ү2/3)) 
RETURN 

PRINT 

GOTO Gl OF 1210,1500,1600 
PRINT "ВОМ №0. 1$"В1 

PRINT 

FIXED 4 

PRINT "COEFFICIENT BUILDUP IS" 
PRINT TAB2"CNU="L1 

BRENT TAB2"CNL="L2 

PRINT TAB10"TOTAL CN="(L1+L2) 
PRINT TAB2"CCF="D1 

BRINT TAB2"CCR="D2 

PRINT TAB10'"TOTAL CC="(D1+D2) 
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1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1415 
1416 
1420 
1430 
1440 
1445 
1450 
1460 
1470 
1480 
1490 
1500 
1502 
1504 
1505 
1506 
1508 
1510 
II 
1595 
1600 
1602 
1604 
1605 
1606 
1608 
1610 
1611 
1612 
1613 
1614 
1615 
1616 
1617 
1618 
1620 
1630 
1640 
1641 
1642 
1643 
1644 
1645 
1650 
1660 


PRINT TAB2''CMU="P1 

PRINT TAB2'"CML="P2 
PRINT TAB2'CMF="P3 

PRINT TAB2''CMR="P4 

PRINT TAB10'"TOTAL CM="(P1+P2+P3+P4)"AT T.E." 
M1= (P1+P2+P3+P4) 

N1=(L1+L2) 

C1=(D1+D2) 

М5-М1-(0.75%М1) 

DISP "INPUT ANGLE OF ATTACK (DEGS.)"; 
INPUT A2 

А2-0.01745%А2 

А9-А2 

N5=N1*COS(A2)-C1*SIN(A2) 
C5sC1*COS(A2)*NI*SIN(A2) 

PRINT 

С2=А2*57.296 

PRINT TAB 2'CL(MEAN)="N5;TAB21''CD(MEAN)="'C5; 
PRINT TAB2"ALPHA ="C2"DEGS" 
STANDARD 

GOTO G3 OF 90,240,250 

END 

N3=L1+L2 

C3=D1+D2 

M3+P1+P2+P3+P4 

FIXED 4 

PRINT '"CN(RMS)*SIN(PHI)='N3 
PRINT "CC(RMS)*SIN(PHI)="C3 
PRINT "CM(RMS)*SIN(PHI)="Mš 
PRINT 

GOTO 250 

N4=L1+L2 

C4=D1+D2 

M4=P1+P2+P3+P4 

FIXED 4 

PRINT ''CN(RMS)*COS(PHI)="N4 
PRINT "CC(RMS)#*COS(PHI)="C4 
PRINT "CM(RMS)*COS(PHI)s'M4 
PRINT 

N1=SQR(N342+N442) 


P1=ATN(N3/N4) 

X=N3 

Y=N4 

Е-Р1 

GOSUB 1900 

Pl=F 

PRINT "СМ(ЕМ5)-"М1І,"РНІ-"Р1%57.296 
C1=SQOR(C3+2+C4+2) 

P2=ATN(C3/C4) 

Х=С3 

Y=C4 

F=P2 

GOSUB 1900 

P2=F 
EM IDEM CCORMS)S"CIS'"PHIsS'P2*57.296 
M1=SQR(M3+2+M4+2) 
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1670 
1671 
1672 
1673 
1674 
1675 
1680 
1690 
1700 
1710 
1720 
15721 
1722 
1723 
1724 
725 
1730 
1740 
1750 
1760 
1761 
762 
1763 
1764 
1765 
1770 
1780 
1790 
1800 
1801 
1802 
1803 
1804 
1805 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1915 
1920 
1930 
1940 


P3=ATN(M3/M4) 

Х-М3 

Y=M4 

F=P3 

GOSUB 1900 

P3=F 

PRINT "CM(RMS)s"M1," PHIs"P3*57.296 
Al1=N1*COS(A9)*SIN(P1)-C1*SIN(A9)*SIN(P2) 
B1=N1*C0S(A9)*COS(P1)-C1*SIN(A9)*COS(P2) 
L9=SQR(A11+2+B142) 

L8=ATN(A1/B1) 

X-A1 

Y=Bl 

F=L8 

GOSUB 1900 

L8=F 
A2=N1*SIN(A9)*SIN(P1)+C1*COS(A9)+SIN(P2) 
В2-М1%С05(Р1)%5ІМ(А9)%С1%С05(А9)%С05(Р2) 
D9=SQR(A2Z+2+B2+2) 

D8=ATN(A2/B2) 

X=A2 

Y=B2 

F=D8 

GOSUB 1900 

D8=F 

A3=M1*SIN(P3)-0.75*N1*SIN(P1) 
B32M1*COS(P3)-0.75*N1*COS(P1) 
М9-50К(А34248342) 

M8=ATN(A3/B3) 

X=A3 
Y=B3 
F=M8 
GOSUB 
M8=F 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT '"CM(C/4)="MS''+"M9"SIN(W1*T+'"M8*57,296")" 
GOTO 90 

IF X«0 AND Y<0 THEN 1920 

IF X>0 AND Y<O THEN 1920 

GOTO 1930 

F=F+PI 

RETURN 

STOP 


1900 


TAB10"AERODYNAMIC COEFFICIENT SUMMARY" 
"CL="NS" +"L9"SIN(W1*T+"L8*57.296"")" 
"CD="C5""+"D9"SIN(W1*T+"D8*57.296")" 
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STEADY DATA OUTPUT 


RUN NO. IS 123 


COEFFICIENT BUILDUP IS 
CNU= 0.3480 
CNL=-0.4281 
TOTAL CN=-0.0802 
СОЕ=-0.0559 
ССК= 0.0721 
TOTAL СС= 050282 
CMU= 0.1092 
CML=-0.2840 
CMF= 0.0017 
CMR= 0.0038 
TOTAL CM=-0.1693 AT T.E. 


CL(MEAN)=-0.0765 CD(MEAN)= 0.0451 (CM)C/4(MEAN)=-0.1092 
ALPHA =-4.9991 DEGS 


UNSTEADY DATA OUTPUT 


CN(RMS)*SIN(PHI)=-1.5385 
CC(RMS)*SIN(PHI)=-0.0312 
CM(RMS)*SIN(PHI)=-0.6336 


CN(RMS)*COS(PHI)=-1.4968 
CC(RMS) *COS(PHI)=-0.0042 
CM(RMS) *COS (PHI) =-0.8581 


CN(RMS)= 2.1465 PHI=Z25. 7/6808 
CC(RMS)= 0.0315 PHI=262.2760 
CM(RMS)= 1.0667 PHI=216.4451 


AERODYNAMIC COEFFICIENT SUMMARY 


0765 + 2.1405 SIN(WI*T+ 225.8305 ) 
эрэ 20451 + 0.1629 SIN(W1*T+ 39.2111 ) 
CM(C/4)=-0.1092 + 0.5836 SIN(W1*T+ 63.0435 ) 
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HPPENDIX E 
PRESSURE CDEFFICIENT DATA WITH STERDY Cu 


RUN ND. Чазй! (и=Й. ИИ 
ANGLE OF ATTACK=-S DEG U=!144 FT/SEC 

PRESSURE А 

THP $ COEFFICIENT ТАР + CUEPETCIENT 
| 1.230 д8 -Й. ИВИ 
2 SE 25 СИЕ 
3 02328 ЗИ -И. ИБЕ 
E ЕНЕ 21 SEIS 
5 Й.ВІБ JA -B.1H7 
Б Й. 7ЙБ 83 zuge 
7 A ARA 3H “ШЕЕ 
Н 2.156 BS тИ ЕБ 
Ч И.ИЙБ 35 a an 
17 тИ 87 “ДЕРДЕ! 
|| -A253 38 pis 
12 us 23 -Aay 
5 de ЧИ "impii 
ІШ ӘЙТЕЙ En ив 
ШЕ ZU BES Ча ТИГІЗ! 
ІБ -Й.БИЧ ШЕ ТИСЕ ЯЗ 
(7) “Йй кыд нм HS 
| В ЭЙЕ | ES SES BH 
15 =Й. ЧЕН ЧБ sura 二 
al URBE a7 A BIE 
2| “ПӘНИ ЫН “чак 
22 ӘЙЕЛІН 19 99 
аз -Й.23И Sa ЗИ Чад 
2 22015 Б иа 
ES =Й АН Va И ЕЯ 
2B sS S ши ales 


Гы 
EST 
l 
t<] 


MAG 
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PRESSURE COEFFICIENT DATA WITH STEADY Cc 


RUN ND. 42322 Сизй.ИНБ5 
ANGLE DF ATTACK=-5 DEB Ц-1ЙЙ FT/SEC 
PRESSURE PRESSURE 
TAP 4 CDEFF | CIENT TAR + COEFFICIENT 
| |. ПАЙ 2B Й.!1Б 
2 2.974 29 -1.171 
E Й. 539 30 -1.й72 
Ч Й.Ч7| 可 -Й. 5ИЧ 
с 0.233 32 И. 378 
Б 2.033 EE -Й. 253 
7 -Й. 185 J4 2.013 
Н -Й.ЧЧВ НЬ 1 UKE 
5 -p.EBB 3E 2.233 
ІЙ -4.818 37 2.26! 
11 -Й.95Б 3B „ИЕ 
|2 -|.И7Ч 39 -0. 105 
13 -|. 215 ug -g. 1g 
1 -1.278 Ч! -Й. ИВЕ 
|5 -|.343 uz -Й. 172 
IE -|.355 ЧЭ -Й.233 
E, -1.372 uu -Й. 22| 
IB -1.317 | US -Й. 258 
19 -|. 236 ЧЕ -Й. 215 
ай -|.18g 47 -Й.278 
aj -3.371 ug -Й. 353 
Ea -5,552 ug -2.177 
23 -2.814 SA -A. J22 
ay -2 HEB 5! 2.118 
25 -1.15Ч Y 7.453 
26 -2.718 33 7.81! 
2 - 2.948 
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PRESSURE COEFFICIENT DATA WITH STEADY Cu 
RUN NO. H2303 Си=И . ЙБЧИ 


ANGLE OF ATTACK=-S DEG i= 100 Pian 
PRESSURE PRESSURE 
TAF + COEFFICIENT THR # 2115 54212 
| |. ИИЙ 28 ы чан 
2 1.387 23 Y. 203 
3 H. 182 ЗИ -H B95 
4 Й 55811 al -Aag 
Б 1323 ай uiam 
B ЙКЕ CE zu. Sl 
1 -И. И | ЭЭ И.ИЙ5 
В 4: 387 ШЕ И.ЙІН 
3 "UBA 36 4.423 
ІЙ Het - SH a 
l| -й.НЇ5 38 -Й.Й31 
2 geel 38 -Й.И97 
13 -ІМИЕК ЧИ = 
IH =| MEA HI -И. 1ЧИ 
[5 и == чад 27:14:13 
ЇБ “АМЕ 515 =Й ЕЙ 
g = [URES Si SH 
ІН cau S pau 
E ES ЧБ Е 
20 SEBS] дй cU SET 
2| а ЕВ HB =Й АН 
Z2 SEIL "3 -0.687 
2 zi HU СИ iras 
24 И ЕВ Si 4.437 
ES EE a 1.117 
25 xcu 3 Й. ТИЯ 
2] иша БУ 
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PRESSURE COEFFICIENT DATA WITH STEADY Co 


RUN NO. 42304 Си=Й .ИЧУВ 
ANGLE OF ATTACK=-5 DEB Ц=1ЙЙ ЕТ/ЧЕС 
PRESSURE PRESSURE 
TAP 4 COEFFICIENT TAP 4 COEFFICIENT 
l 1.872 2B 1.237 
2 1.005 29 1.192 
3 1.903 ЗИ -Й.238 
Ч J AE З! -2. 287 
5 1.538 ча -Й. 136 
E 1.312 33 -2. 164 
7 "Meis 34 -Й.ЙІВ 
Н -Й. 197 BE -Й.ИЗБ 
9 -Й. 352 3E -й.И5Ч 
I] -9.519 E -Й. 175 
l| -Й.БІЗ 3B -2. 125 
[2 -Й.775 38 -Й. 119 
13 -Й.928 ЧИ -йЙ.1Ч3 
1Ч -1.01Й TT -й. 167 
E -|.И92 42 -].z8u 
IE - | .B2E 43 -Й. ЗБ 
[7 -1.И5Ч yy -Й.3Б7 
IB -|.Й|Н ug -日 .367 
19 -Й.ВБ7 ЧЕ -Й.3Б7 
22 -7.953 47 -0.534 
21 -|.949 ug -2.777 
22 -2.994 ug -Й.БІЗ 
23 -1.534 СИ -Й.297 
au -Й.ИЧЕ 5! -0.379 
25 -日 .229 G2 Sy 
25 -Й.782 E 1.579 
27 -Й.ИИВ 


ЕЕ UU ORE EE An tst 





= 
= 
“г 
de 


ANGLE OF ATTACK=-S DEG 


ESSE 
PEREAS TENT 


а ШЕ шеш — 


nd ЙИК Л ПГ AAA = =— = 


ERE URE CME ЕМ DATA NETA STERDY <u 
RUN NO. 42345 


ЕН 
. HBH 
983 
.HIB 
Ez 
ae 
. 268 
. ZB7 
as | 
are 
‚АИЯ 
‚ БИЧ 
RE 
887 
ЗЕ 
ТЕБЕ 
. BBB 
AS 
Bal 
ДЕЛ) 
E 
2 
ЕНЕ 
ДЕС 
221415 
‚ ЙЗИ 
MBE 
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ІНЕ % 


28 
2H 
ЗИ 
ا3‎ 
32 
33 
34 
38 
36 
37 
38 
39 
ЧИ 
Ч | 
42 
43 
HH 
HS 
HE 
47 
HB 
-3 
ЕЙ 
x 
52 
53 


Cu= W272 


(= ПАНЕВ ЗЕ 
PRESSURE 


НЕРІ ЛЕНІ 


E 
BS 
‚ИБИ 
‚ИЧЕ 
‚232 
‚ИВР 
‚ИЧА 
MES 
¿287 
. 128 
.|BE 
‚224 
224 
‚22y 
250 
396 
‚432 
.HS | 
44g 
629 
‚898 
-ВЧ2 
‚749 
‚533 
‚262 
.183 





PRESSURE ОШЕН БЕ С ЕЧЕНТН НЇН STEHOY C. 


RUN ND. 42328 C,-240.8133 
ANGLE OF ATTACK=-5 DEG И=ТИЙ FT/SET 

PRESSURE PRESSURE 

ТЕР # “ШЕЕ ВМВ THP $ CHEE CENT 
| a S22 2 日 Й.ИНН 
2 4.2948 24 2.285 
3 Й 5H ЗИ 1.275 
4 Й.НДН z 4.463 
5 1. 743 BZ 4.439 
Б 1.534 аа “Ae dae 
D 4.318 E -И. ИВИ 
В Й. ИЧЕ J5 -2. 152 
= В ЗБ -И.22й 
ІЙ -ШЕЗ | | Ej “4.273 
|| ӘЛІН? 3E аа 
ТЕ ЗІНЕ Ба zy 25 
ES -Й.БИЧ ЧИ -4. 25! 
а -0P E97 51) =Й че 
|5 zur 789 2 zu 
ЇЕ ABE нэ “=Й ЫЕ 
17 =A BSH КІН US 
IB "UNE a zu 363 
3 mud HE AB 
zu SA BH! ail -5.772 
21] A 375 HB e s 
dà E a9 IES 
аз = ЧИЕ KH zu ЕН 


24 Й.ЙЙБ Sl 21201205 
AS 4.467 SZ -4.B3B 
26 4.297 E Bras 
27 2.267 
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ТЕР # 


O J Ш ий 2 U — 


N hi N J pJ — — — — — — — — — — 
cu — El EH X im DICE mi Ar — ыш 


АБ 
2E 
z 


“сг 
. Las 
‚ЧИН 
1812 
ЛЕЙ 
БЕН 
- 389 
El 
. ИЧВ 
ZEE 
„274 
2 
. ЧВИ 
Sa 
.E2E 
Жа! 
¿535 
S87 
. 138 
Hall 
987 
li 
eis 
208 
25 
. 33H 
НАҚ 


APPENDIX < 
UNSTERDY PRESSURE DATA 
ANGLE OF ATTACK=-S DEG 


MEAN PRESSURE 
COEFFICIENT 
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Ц=1ЙИ FT/SEC 


UNSTERDY PRESSURE 
COEFFICIENT 


га 


м мы NN NH NS NANA NH NAN NN NA NH ANN ы 


-HAH 
. 186 
HIS 
. БВ 
. 203 
2211 
NEZA 
. 38H 
. ЧИЙ 
SIS 
Ala 
. 68 | 
Ma 
ME 
3513 
. ]BB 


ES 


.ББІ 
. БИЙ 
SEZ 
ЕЕЕ! 
.ЧИЙ 
. 338 
TABS 
args 
. 188 
.| 日 日 


PHASE 
HNBLE 
| БИ 
| НИ 
ВУЙ 
Ваи 
ӘНЕ 
ERES 
И 
| 7 
23 
аа 
ач 
37 
HB 


ug 
= 


Б | 


SS 
ES 
БИ 
БН 
75 
BB 
| ИЗ 
lll 
| ИН 
ум 
ЫН 
BR 





THP 2 


2H 
29 
ЗИ 
EI 
a2 
33 
En 
EN 
EE 
37 
3H 
E 
ЧИ 
HI 
2 
юэ 
ki 
55 
ЧЕ 
gl 
ЫН 
HH 
EVA 
EM 
ne 


UNSTEADY PRESSURE DATA 
C CONTINUED 9 


ANGLE DF ATTACK=-S DEB U=144 FT/SEC 
MEHN PRESSURE UNSTERDY PRESSURE PHASE 
CDEFFICIENT COEFFICIENT HNBLE 
-Й.ЧІЙ Йе ды 348 
-j . 285 4.434 ASA 
i528 И. |Б] АНН 
UNES A E23 74217 
-H.1BI 14.311 292 
СИРИЙН 4.376 zus 
24 178 ا‎ 2724 
de sie И. Ч7В 255 
БИ УЗИБ 4.539 zyj 
4 22H И. БИИ 247 
Seeds j. SBE 232 
-И.26Й И.БВН ва! 
A 28% И.БНЧ 2040) 
SH И.БНИ 224 
EE ӘСІ Zum 
-j . SHE шз ZIS 
SA SSH И.БИЙ alí 
=Й aos 1.538 АЙК 
-И.БЧН 4.513 ISI 
4:853 4.963 [ЧИ 
= ЛЕЙ MIES ІНЕ 
Е ДЕРДЕ |75 
zA SAJ 4.445 17! 
= An Sl [Sa 
1125 И. СЯН 157 
=й. чо zs 12 


Sg 


— s мш 
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